Low-frequency dynamics of disordered XY spin chains and pinned density waves: 
from localized spin waves to soliton tunneling 
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A long-standing problem of the low-energy dynamics of a disordered XY spin chain is re-examined. 
The case of a rigid chain is studied where the quantum effects can be treated quasiclassically. 
It is shown that as the frequency decreases, the relevant excitations change from localized spin 
waves to two-level systems to soliton-antisoliton pairs. The linear-response correlation functions are 
calculated. The results apply to other periodic glassy systems such as pinned density waves, planar 
vortex lattices, stripes, and disordered Luttinger liquids. 



In this Letter I revisit the problem of the low-frequency 
response of a periodic elastic string pinned by a weak ran- 
dom potential. This problem has been repeatedly studied 
in connection with dynamics of density waves, Luttinger 
liquids, glasses, and random-field XY spin chains. It 
is also germane for statistical mechanics of periodic pla- 
nar systems (Josephson junctions, stripes) with columnar 
disorder. The approach presented below is customized 
to one dimension (ID) but the results should apply to 
all dimensions d < 4 after a few minor modifications. 
The difficulty of the problem stems from strong nonper- 
turbative disorder effects at low enough frequencies and 
large enough length scales. Such effects are termed pin- 
ning [1] or localization depending on a context. Rigor- 
ous results regarding the pinned regime are restricted to 
one particular value of the chain elastic modulus where 
the problem maps onto noninteracting fermions in a ran- 
dom potential, whose dynamics is governed by the fa- 
mous Mott-Halperin-Berezinskii (MHB) law [2, 3]. In 
this special case the chain is extremely soft and behaves 
as a quantum object. On the other hand, real charge 
density waves (CDW) and low-density electron gases are 
much more rigid and consequently quasiclassical. Despite 
many previous attempts [4-6] to calculate the response 
of the latter type of systems, no fully satisfactory so- 
lution has emerged. The existing controversy between 
different authors is clouded by an uncontrolled nature of 
approximations they use. The purpose of this Letter is 
to demonstrate a feasibility of a quantitative analysis of 
the problem and to elucidate the physical picture of the 
low-frequency excitations in the glassy phase of a rigid 
(strongly interacting) system. Most of the discussion is 
focused on the classical dynamics; however, at the end I 
also consider quantum effects and establish the connec- 
tion with the MHB law. In plain words, I offer a definitive 
answer the following basic question: // we shake a pinned 
elastic chain, how does it respond? 



Model. — Consider an XY spin chain with spin stiff- 
ness c, spin wave velocity v, lattice constant a, which is 



described by the Lagrangean 



(1) 



where Xj = aj and hj are quenched Gaussian random 
variables with zero mean and variance = A. De- 

pending on the context, can equally well represent the 
phase of a CDW, the bosonized density field of a Lut- 
tinger liquid or a transverse displacement of a flexible 
line object. The quantity of interest is the propagator 
D{uj, q) of the field or equivalently, the quantity 



ct(w,(?) = -i 



47r2n2 



(2) 



which has the meaning of the ac conductivity. Of partic- 
ular interest is a{uj,q = 0) = cr(w). 

Classical glass. — A classical chain respond to an exter- 
nal perturbation by elastic vibrations around the ground 
state 0°(x). It is therefore instructive to recall the basic 
properties of 0°, described by the collective pinning the- 
ory [1, 4]. The spatial structure of 0° is determined by 
the competition between disorder and elasticity: on short 
scales the elasticity prevents large phase distortions but 
on long scales the disorder eventually wins and the dis- 
tortions grow without a bound: ([0°(a;) — 0"(O)]2) oo 
a.s X oo. There is a characteristic scale Rc where a 
27r-distortion is accumulated. At this scale the typical 
elastic and disorder energies of a chain segment are of 
the same order, c/Rc ^ \/ RcJ^/a. It is easy to see that 
Rc ^ (c^a/A)^/'^. A crude but useful picture [4] is to 
imagine that the chain consists of domains of size Rc in- 
dividually pinned by collective potential wells generated 
by random fields hj . 

An explicit algorithm for finding 4P{x) was originally 
given by Feigelman [5] . Consider the segment of the chain 
with j leftmost spins, which satisfies the boundary condi- 
tions at the left and has the angular coordinate of the jth 
spin fixed at a given value (j)j . Let [xj , 4>j ) be the min- 
imal energy of such a segment optimized with respect to 
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all spin coordinates in the interior of the segment. Func- 
tion E~ satisfies the recurrence relation [5] 

E-{x,+i,q^) = min{E-ix,,q^') + ^i^-^'f} 

+ aU{xj+i,(l)), (3) 
?7(x,+i» EE a-\h,+,e^^ + h*^,e-^^). (4) 

Let E^ {xj ,(f>) be a similar function for the right end of 
the chain, then the desired (f>^{xj) is the value of that 
minimizes the sum E(xj , (jj) = E~ (xj , (f>)+E~^ (xj , 4>) . For 
Rc ^ a, (fP varies little from j to j + 1, and so Eq. (3) 
possesses a meaningful continuum limit, 

E- =-{ll2c){E-f + U{x,4>), (5) 

which is the Kardar-Parisi-Zhang (KPZ) equation [5, 7]. 

Due to the 27r-periodicity of U the solutions of Eqs. (3) 
and (5) become 27r-periodic at large x irrespectively of 
the boundary conditions at the ends. A typical behav- 
ior of E{(j)) is illustrated in Fig. 1. £'(0) gives direct 
information about the rigidity of the system. The min- 
imal work needed to twist a given spin of the chain to 
the angular coordinate is equal to £'(0) — E{ip'^). For 
small 6(f> = — 0°, it is quadratic in 50. The elastic 
distortion caused by the twist is localized predominantly 
within a domain of length Rc around the chosen spin. 
When |50| exceeds a certain value, the chain suddenly 
snaps into a conformation corresponding to a compet- 
ing metastable state. At such 0, function i?(0) exhibits 
upward cusps (typically, one per period — see Fig. 1), 
extensively studied in the context of the KPZ equation 
and Burgers turbulence [7, 8]. 

Of interest to us here are the low-frequency excitations 
of the chain. They can be visualized as localized mechan- 
ical oscillations of essentially rigid segments of size Rc- 
This concept can be succinctly expressed by means of the 
following low-frequency local effective action 

- ^MiMf - £(0), (6) 

which describes the motion of a particle ("domain") of 
mass M ~ Rcc/v^ in the potential E{(j)). To justify this 
action, I will follow Ref. [5] but with important modifi- 
cations, leading to very different end results. 

To calculate the linear response we expand £ [Eq. (1)] 
in (50 and keep only the second order terms. After that 
we perform the usual spectral decomposition to obtain 

D(u;;x,x ) ^ - y -—^ — , eiu)) = ^r-. (7) 

^ ^ a ^ En - e(Lu) + lO ^ ^ 1)2 ^ ^ 

The frequency aj„ of nth eigenmode and its wavefunction 
■^n satisfy the discrete Schroedinger equation 

-|V^l/'n(a;j) -f lJ{Xj,(\P(Xj))-^ri(Xj) = £„-0n: (8) 
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FIG. 1: Typical behavior of E{(fy (for a given fixed x). 

where V is the lattice derivative and e„ ~ siyJn) is the 
"energy" of the mode. From the scaling of the "kinetic" 
and the disorder energies with distance, we quickly de- 
duce that the low-energy modes are necessarily trapped 
in the potential wells of the random potential \] . To 
determine their typical localization length Lioc we can 
compare the kinetic energy c/L"^^^ of a wavepacket of 
size L\oc with the depth A/LiocO of a typical poten- 
tial well, which yields iioc ~ Rc- Thus, Rc is the 
unique characteristic length of the classical glass regime, 
which shows up both in statics and dynamics [9]. The 
crossover to this regime occurs around the pinning fre- 
quency LOp — vl Rc- Note that the depths of individual 
potential wells have a broad distribution around the typ- 
ical value A/ Rett ~ e(i^p), which gives rise to an in- 
homogeneously broadened spectrum extending down to 
£ = 0. Modes with e„ ^ s{u!p) stem from near cancella- 
tions between the kinetic and potential terms. 

Let X be a point near which one of such modes, ipn, 
is localized, then for all ut < ujp, Eq. (7) implies that 
D(u!; X, x) = A/ (e„ — e + iO) + smaller terms, with A = 
V4(a;)/a Rc^- On the other hand, [D{0;x,x)]-^ = 
E^^{(jP) = a{x) [5], and so [D{lo] x, x)]~^ = a- Mu)^ + 
o{a)- This indicates that CeS indeed has the correct form 
for small (50. (This is all we need to start using £eff for 
calculating the classical linear response). Equation (6) 
is also correct for adiabatically slow motions, dt(j) 0. 
Thus, it cannot be badly wrong for all (50 smaller than 
the distance from 0o to the nearest cusp of i?(0). In- 
deed, local smooth distortions should appear to the rest 
of the system as adiabatically slow ones provided their 
characteristic frequencies are sufficiently small, u <^ Up. 

From Eq. (6) we see that the disorder-averaged low-w 
response of the chain is encoded in the small-(a; behavior 
of the distribution function P{a) of the local ground- 
state rigidity a(x), e.g., the spin wave density of states 
p{e) = {Nsa)~^ J2n - ^ri) is given by 

p{e) - P{a = eRc), e < e(wp). (9) 

Here Ng is the number of spins in the chain. Similarly, 
starting from the spectral representation 

'^^^i^) = ^{T.ndlSie-e^)), (10) 
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FIG. 2: Symbols: p{e) and P{a) computed and averaged over 
10® disorder realizations numerically. Siirmlation parameters: 
c=l, Ns = 200, A = 7 X 10"'^. Taking previous definitions 
literally, we get Rc ~ 112a. The actual correlation length of 
^° is about 10a. Straight lines: fits to the power laws with 
exponent s = 1.7. 

where d„ = ^jtpn{xj), we obtain d„ <^ {Rc/aY^^ and 

3ica{uj) ^ e'^ujRcf)[e{uj)], uj < ujp. (11) 

Now I intend to show that at small a 

P{a) (X a', (12) 

where s = 3/2. Combined with Eqs. (9) and (11), this 
entails [10] 

^ea{uj) ~ {e'^vRc/c){ijj/ujpY, u) <C ujp. (13) 

Recall that a is defined to be the second derivative 
of E[(l)) at the point of its global minimum (/)o. Let 
us first demonstrate that at arbitrary local minima 

Prob(Q!) = Pi{a) oa \a\ at small a. Indeed, lo- 
cal extrema are the points where E^{(j)) ~ 0; hence, 
Pi{a) = (AT-i d4)5{E^^-a)5{E^)\E^^\) where iV„ = 

d(()5{E^)\E^^\. (The term \E^^\ in the integrands is 
the Jacobian). In the limit a ^ we obtain 

Pi{a) = \a\ (n-^ Jd4>5{E^^)5{E^)^ = Co|a|. (14) 

The main contribution to the integral that determines 
the constant Co ^ R^/c^ is supplied by configurations 
where i?^^ and E'^^ have their typical values of the order 
of c/Rc but are opposite in sign and almost cancel each 
other. Thus, soft modes arise from the frustrations in 
the ground state, e.g., when cj) = (ji^ gives a rather low 
energy to the left half of the chain but corresponds to a 
high energy state of the right half. Such frustrations can 
always happen because E~ and S+ are uncorrelated. 
Let us now consider the global minimum. In its vicinity 

E{^) = E{<}P) + ^5<j>'' + f3Sct>^ + ^5,i>\ (15) 



which leads to the double-minimum structure depicted 
in Fig. 1. It is easy to see that is the lower of the 
two minima only if |/3| < ^flarj. Via a straightforward 
extension of the argument leading to Eq. (14), one can 
obtain Pj(a,/3) = Ci\a\ for the joint distribution func- 
tion of small a and /3. The restriction on allowed /3 at the 
global minimum leads to P{a) cx ^/aPi{a, 0) cx a^/"^ as 
I claimed above. Here I ignore competing minima away 
from the immediate vicinity of cjp because their number 
within a 27r-period is Nm ^ 1. They cannot yield any 
additional powers of the small parameter a. 

Now I present my numerical results. The numerical 
solution of Eq. (3) does not pose major difficulties. This 
has to be done for a number of disorder realizations {/ij}, 
followed by solving the eigenvalue problem (8), calculat- 
ing the desired response functions, and their statistical 
averaging. Such a procedure produces the rigidity distri- 
bution P{a) shown in Fig. 2. The power-law behavior 
predicted by Eq. (12) is apparent. Function p(e), plot- 
ted in the same figure, can be well fitted to the power 
law with the same exponent, which shows that Eq. (9) is 
also well satisfied. The actual value of the exponent from 
those fits, s = 1.7, deviates slightly from my analytical 
prediction s = 3/2. Further numerical work is needed to 
resolve this small discrepancy but there is enough accu- 
racy to rule out s = 1/2 [4, 6] or s = 1 [5]. As for the 
functional shape of E{4>), I observed the double- minimum 
structure sketched in Fig. 1 quite often, roughly in 20% 
of the cases. 

Two-level systems. — For a typical double- well poten- 
tial E{<j)) with a given oscillator frequency uj, the dis- 
tance in (j) and the energy barrier between the two min- 
ima are of the order of uj/u)p and {c/Rc){uj/ujpY, respec- 
tively. Both decrease with uj, and it can be verified that 
at a; ~ = K^^^ujp where K = hv/2iTC <C 1, the matrix 
element / for the quantum tunneling between the two 
minima becomes of the order of ftw. This implies that at 
frequencies below ujq the response of the chain is domi- 
nated by quantum tunneling. It contributes to ^ea{w) 
whenever the levels localized in the two minima are split 
in energy by exactly hw. A straightforward analysis [11] 
of the model defined by Eqs. (6), (14) and (15) leads to 

2 

^ea{uj) ~ -tRcK'^ — , Wgi « w < ojq. (16) 
ii (jjp 

The origin of the frequency scale Uqi = ujpexp{—K~^) 
is as follows. The energy splitting of our two- level sys- 
tems (2LS) is bounded from below by /. Therefore, as 
(jj decreases, the tunneling barrier and the tunneling dis- 
tance have to increase. At w < LVgi a typical 2LS cannot 
possess such a small I and Eq. (16) ceases to be valid. 
As w continues to decrease, the dissipation is initially de- 
termined by rare 2LS with unusually large barriers, until 
at LUg = LUp cxp{—K~^/^) the soliton mechanism becomes 
more prominent. 
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FIG. 3; The sketch of ^ea{uj). The nonperturbative low- 
frequency regime uj < Up is discussed in the text. The formula 
for Lu > Up where the perturbation theory applies is from [13]. 



Soliton tunneling. — In the soliton mechanism a large 
tunneling action (small /) is due to a large tunneling 
mass. The object that tunnels is a segment of the chain 
of length ^tun(^) 3> Rc- The optimal way to accom- 
plish such a tunneling is to send a virtual soliton (a 
27r-kink in (j)) over the distance Ztun. Another route to 
the concept of solitons is appealing to universality of the 
low-frequency physics of the glass phase (throughout the 
range [13] K < 3/2). If it holds, the conductivity must 
be calculable by generalizing the Mott-Halperin argu- 
ment [2]. This argument devised for = 1 focuses on 
the tunneling of single electrons. But they are precisely 
the 27r-solitons in the bosonized formulation (1), so re- 
placing the word "electron" by "soliton" everywhere in 
the argument seems entirely natural. 

The ID MHB law can be written in the form [14] 

^ea{Lu)r^Q'lljspUhw), (17) 

where Ps{E) is the density of states, /tun ~ ln(/o/Sii') 
is the typical tunneling distance, and Is is the localiza- 
tion length of the tunneling charge-Q objects {Q = e for 
27r-solitons). A typical soliton has the size Rc and the 
creation energy Eg = c/ Rc, which implies Ig ~ hv/Eg ~ 
KRc [15]. Under the standard assumption [16] that 
Ps{E) ~* const ^ (RcEs)^^ for i? ^ 0, we get 

2 2 

^ea{Lu) ^ ^RcK-''^lii^ uj-^Us. (18) 
h uj^ ujp 

The frequency dependence of the conductivity is sum- 
marized in Fig. 3. My results support the notion that 
the infrared behavior of nondissipative classical [11, 12] 
and quantum [13] glasses are universally controlled by the 
w^-dependence (13) and the MHB law (18), respectively. 
In the model considered, as K increases, the system be- 
comes more quantum. The quasiclassical regimes (2LS 
and oj'^) shrink and eventually get eliminated at if ~ 1, 
where the soliton tunneling crosses over directly to the 
Drude behavior, in agreement with Ref. [3]. 



Relation to experiments. — The obtained results have 
implications for a wide variety of physical systems. Per- 
haps, the most studied of them are the CDWs. Evi- 
dence for the 2LS in CDW compounds has been seen in 
the specific-heat [17]. Transport data typically fit the 
a{uj) (X Tu!^ dependence with 6 « 1 [18]. It may or may 
not be due to the 2LS physics, but in any case, it means 
the true zero-temperature limit I studied here has not 
been realized in the experiment. To verify Eqs. (13) and 
(16) new ac transport experiments at considerably lower 
T are needed. As for Eq. (18), strong modifications due 
to Coulomb effects [14] are expected in real CDWs. An- 
other class of materials to test would be the high-spin 
spin-chain compounds whose dynamics can be studied 
by electron spin resonance. Finally, some of the ideas 
behind Eq. (13) may also be relevant for understanding 
the w^-broadening of phonons in glassy liquids [19]. 
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